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1 Introduction 


The propagation of an evanescent wave through a potential barrier has long been studied 
[H121 El St The question of tunneling time i.e. the time spent by a wave packet with mean 
incident energy smaller than barrier height in the classically forbidden region and the anal¬ 
ogous situation of evanescent waves in optics have in recent years attracted considerable 
attention. One of the most interesting aspect of these studies is the saturation of tunnel¬ 
ing time with respect to the width of the barrier and is referred as Hartman-Flctcher (HF) 
effect. Hartman studied [5] the tunneling time by constructing a metal-insulator-metal 
sandwich by the method of stationary phase to demonstrate the experimental agreement 
of the in dependency of tunneling time on the width of the barrier. Fletcher indepen¬ 
dently showed [6j the saturation of time delay by considering tunneling of evanescent wave 
through a thick barrier. This exciting result triggers series of famous experiments with 
incident wave both in microwave m-0 and optical range HDl-HU to obtain the saturation 
of tunneling time by considering single, double, and multiple real barrier potentials. Re¬ 
cently Longhi et. al HU have measured tunneling time for a double barrier optical grating 
and found that the tunneling time is paradoxically short and independent of barrier width 
and separation between the barriers. Similar conclusion was obtained by Dutta-Roy et. 
al (12] by considering single barrier associated with inelasticity. 

In the present work we study the HF effect by considering an array of non-Hermitian 
barrier potentials. The motivation arises from the huge applicability of non-Hermitian 
system [13]- [15] in various branches of physics over the past one and half decades. In 
particular non-Hermitian theories are considered to be the topic of frontier research work 
on transport [16j-[18| and scattering [I2]-|3Uj phenomena for matter as well as electromag¬ 
netic waves. Scattering from complex potential has very rich features and leads to many 
technological developments [26 ] I3T] 32] , Even though the HF effect have been extensively 
studied both theoretically [33j-[36j and experimentally 0-n for real periodic potentials, 
it has not been discussed for complex potentials which are associated with both elastic 
and inelastic channels, except the work in [Q] where the tunneling time for a single com¬ 
plex barrier has been calculated and HF effect is discussed for the weak absorption. For 
the strong absorption the HF effect for a single complex barrier is shown to disappear. 
This result is consistent with the experimental findings in Ref. [7]. Being a quantum me¬ 
chanical process the characteristics of tunneling do not guarantee to be the same when 
wave transmits through multiple number of barriers which are separated from one an¬ 
other. Therefore it is worth investigating the various characteristics of tunneling time 
and the presence of HF effect when waves scatter through an array of complex barriers. 
We consider an array of non-Hermitian square barrier potentials with a fixed height and 
periodicity which is actually the mathematical structure of the famous Cronning-Penny 
model that depicts the lattice structures in crystal. In this work we calculate the time 
taken by a wave packet to tunnel through such an array of complex barrier potentials to 
show the various characteristics of HF effect. We put some light on the the behavior of 
tunneling time and discuss the long debated HF effect i.e. the saturation of tunneling time 
for matter waves scatter through such arrays of non-Hermitian potentials. We calculate 
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the tunneling time by using the method of stationary phase for an array consist of any 
number of complex barrier potentials and discuss the saturation of tunneling time with 
the number of the complex barriers by varying different parameters in the system. The 
saturation crucially depends on the coupling V c between elastic and inelastic channels. 
Saturation is achieved only for small V c . This implies system shows HF effect only for less 
absorption which is consistent with the results in mm ■ For fixed V c HF effect depends 
on the width of the individual barrier and separation between the consecutive barriers. 
Unlike the array of real barriers here the saturation is obtained only for certain range of 
width of individual barrier. One of the rich features of tunneling time for real barriers 
is the occurrence of resonances [37] for specific values of the separation between adjacent 
barriers. We observe some of resonances disappear with increase of absorptivity in the 
system of array of complex barriers. We show that the resonances in tunneling time for 
array of complex barriers is regulated with increasing V c . We further demonstrate the 
existence of HF effect even for more realistic situations where inelasticity is randomly 
chosen for individual barriers. In case of emissive inelastic channel we show that the HF 
effect occurs only at certain small values of incident energy. 

Now we present the plan of this paper. We review the HF effect for single barrier in 
Sec.2 to outline the methodology for later sections. In Sec.3 tunneling time is calculated 
for an array of complex barriers. Various results regarding the HF effect for arrays of 
complex barriers are discussed in Sec.4. Sec.5 is kept for conclusion and remarks on 
further prospectives. 


2 HF effect for single barrier 

In this section we review the method of calculating the tunneling time of a wave packet 
through a single real barrier by stationary phase method [38]. In this method the tunneling 
time (r) is defined as the time taken by the peak of the incident wave packet to traverse 
the classically forbidden region and emerges as transmitted wave packet. To calculate r 
we consider the evolution of an incident localized wave packet which is described by, 

J Gko(k)exp(ikx — iEt/h)dk . (1) 

Gk 0 (k ) is the normalized Gaussian function of wave number k peaked about the mean 
momentum hk 0 . Due to presence of a barrier an incident wave packet would emerge after 
transmission as 

J Gk 0 | clt | exp(ikx — iEt/H + i$(k))dk (2) 

where clt is the transmitted amplitude ax —\ dr(k) \ e z ®( k \ The time r at which the peak 
of the wave packet emerges from the barrier (V = V r for 0 < x < b and V = 0 elsewhere) 
is the given by stationary phase method as, 

-rr[ikb — iEr/h + i$(k)] — 0 . (3) 

CxjKj 
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Figure 1: Tunneling time with the width of the barrier shows Hartman-Fletcher effect for 
a real barrier with V r — 1, E — 0.5 and h — 1 — 2m. 


This implies 


r = h —r^ + 


dE ' hk/m 

For a single real barrier the transmission amplitude is given by, 

2 kqe~ ibk 


(4) 


ar(k) = 


2kq cosh(bq) — i(2E — V r ) sinh(6g) 
Thus r is calculated in stationary phase method 


with k = \p2nxE, q = \j2m{V r — E)/h (5) 


as, 


„ d 
t = n—— tan 
dE 


-l 


k 2 - q 2 
2 qk 


tanh qb 


( 6 ) 


From the above equation we observe that r —> 0 as b —> 0 as expected, however when 
b —y oo, r = i.e. tunneling time is independent of the width of the barrier b for 
sufficiently opaque barrier. This paradoxical result (HF Effect) is demonstrated in Fig. 

m 

The Fig. [1] shows the saturation of tunneling time when width of the barrier is larger 
than a certain value and incident energy is less enough than the height of the barrier. 
Recently the tunneling time has also been calculated for single complex barrier potential 
to show the HF effect JT2J for weak absorption. HF effect disappears when the imaginary 
part of the barrier has large value. In the later sections we discuss the tunneling time and 
HF effect for array of complex barriers. 
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3 Tunneling time for array of complex barriers 


In this section we intend to calculate tunneling time for an array of complex barriers. For 
this purpose we start with a single complex barrier with a non-Hermitian potential, 

V(x) = V r + iVi for ai < x < d 2 

= 0 elsewhere , (7) 


where V r and Vi are real and a 2 — a\ = b is the width of the barrier. Following the idea of 
[12, 39j this complex barrier can physically be realized by two channel formalism, where an 
absorptive (or emissive) inelastic channel is evanescently coupled with the elastic channel 
via a coupling potential V c . The real part of the potential is corresponding to an elastic 
channel described by the Schroedinger equation (with h — 1 — 2m), 


d?_ 

dx 2 


+ V r 


-0( x ) + V c (j){x) = Ei/j(x), 


( 8 ) 


Whereas the imaginary part of the potential is associated with an inelastic channel and 
is described by the Schroedinger equation as, 


dx 2 




<f>{x) + V4>(x) 


(.E — A)ijj(x) 


(9) 


A is the energy absorbed by the system from incident wave due to the presence of the 
imaginary part of the potential. For positive A, the inelastic channel is absorptive whereas 
it is emissive for negative A. The scattering wave functions in the elastic and inelastic 
channels (as in Ref. 112) can be obtained by solving Eq. (JHJ) and (J9]) successively, 

= ( Be ax + Ce~ ax ) sin(0/2) + {Fef* x + Ge~ Px ) cos(0/2) 

<j>{x) = {Be ax + Ce~ ax ) cos(0/2) - {Fe Px + Ge~ px ) sin(0/2) (10) 

where, 

a 2 = V r — E + ——^—^-(secd — 1) ; /3 2 = V r — E — ——^—^-(sec6* — 1) ; (11) 

md 6 = tan " U- 2 k- a ) (l2) 

The asymptotic forms of these scattering wave functions in the elastic and inelastic chan¬ 
nels are given as follows, 

at x < ai, ^~{x) = Ae lkx + De~ lkx ; 4>~{x) = R ine ie~ lkx , 

at x>a 2 , ^ + {x)=Pe ikx + Qe~ ikx ; <j) + (x) = T inel e ik ' x (13) 

where k = \f~E, k 1 = \JE — A. As the wave functions of both the channels in Eq. (flOl) , 
and their derivative satisfy the condition of continuity at x > a± and x < a\, one get total 
8 equation for the two individual channels which are, 

Ae iaik + De~ iaik = (Be aai + Ce ~ aai ) sin(0/2) + (Fe^ ai + Ge~ pai ) cos(0/2) (14) 
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ik{Ae iaik - De~ iaik ) = a (Be aai - Ce ~ aai ) sin(0/2) + /3 ( Fe 0ai - Ge~ 0ai ) cos(0/2) (15) 
Pe ia2k + Qe~ ia2k = (Be aa2 + Ce ~ aa2 ) sin(0/2) + (Fe 0a2 + Ge ~ /3 “ 2 ) cos(0/2) (16) 

ik(Pe ia2k + Qe~ ia2k ) = a (Be aa2 - Ge" Qa2 ) sin(0/2) + p ( Fe 0a2 - Ge~ 0a2 ) cos(0/2) (17) 


Rineie~ iaik ' = (Be aai + Ce~ aai ) cos(0/2) - ( Fe 0ai + Ge~ 0ai ) sin(0/2) (18) 

- ik'R inel e~ iaik ' = a (Be aai - Ce~ aai ) cos(0/2) - p (Fe flai - Ge"^ 1 ) sin(0/2) (19) 

T ine ie ia2k ' = (Be aa2 + Ge" 002 ) cos(0/2) - (Fe^ 2 + Ge~ fiaa ) sin(0/2) (20) 

ik'T inel e ia2k ' = a (Fe aa2 - Ge" aa2 ) cos(0/2) - p (Fe^ 2 - Ge“** 2 ) sin(0/2) (21) 

During a bidirectional scattering in the elastic channel the asymptotic amplitudes in 
the right hand side (i.e. F, Q) and those for the left hand side (i.e. A, D ) are related to 
each others via the transfer matrix as, 


( P \ = ( M u M 12 \(A\ 

[q ) [m 21 m 22 ){d ) 


( 22 ) 


From Eqs. (TT41) - (f2Tj) and by using Eq. (1221) we calculate all the components of the M- 
matrix as, 


M u = 


0 i(ai — a 2 )k 


Akafi 


a{ik 2 z + 2 ky/3 — iz/3 2 ) cos 2 (6 l /2) + P{ik 2 v + 2 kwa — iva 2 ) sin 2 (0/2) 


i{kk'za — ( kk'v + i(k — k'){w — y)a + va 2 )f3 + zap 2 ) 2 sin 2 6 
A(3(k' 2 v + 2 ik'wa — va 2 ) cos 2 (0/2) + Aa{k' 2 z + 2 ik'y/3 — z/3 2 ) sin 2 (0/2) 


(23) 


M12 — 


g— i(ai+a2)fc 

Akot/3 


—iza{k 2 + /S 2 ) cos 2 (0/2) — iv{3(k 2 + a 2 ) sin 2 (6 l /2) 


+ 


i (fc 2 {F(^q( — u/3) — i(u) — r/)a/3} 2 — a 2 p 2 {va — z/3 — ik'(w — y)} 2 ) sin 2 6 
A/3(k ,2 v + 2 ik'wa — va 2 ) cos 2 (6>/2) + Aa{k' 2 z + 2 ik'y/3 — z/3 2 ) sin 2 (0/2) 


(24) 


M 2 i — 


s i(ai+a 2 )fc 


4/ca/3 


iza[k 2 + /3 2 ) cos 2 (0/2) + iv/3(k 2 + a 2 ) sin 2 (0/2) 


— i (jv 2 {k'(za — v/3) + i(iu — ^cqS} 2 — cx 2 /3 2 {ra — z/3 — ik'(w — y)} 2 ) sin 2 9 
A(3(k ,2 v + 2 ik'wa — va 2 ) cos 2 (6 l /2) + 4a(F 2 £ + 2 ik'y/3 — z/3 2 ) sin 2 (6 l /2) 


(25) 


M 22 — 


e i(a 2 -a 1 )k 

Aka/3 


a{—ik 2 z + 2 ky/3 + iz/3 2 ) cos 2 (6*/2) + (3{—ik 2 v + 2/cu>a + iva 2 ) sin 2 (6*/2) 


iikk'za + {—kk'v — i{k + /c')(w — y)a + va 2 )(3 — za/3 2 ) 2 sin 2 6 
A/3(k' 2 v + 2 ik'wa — va 2 ) cos 2 (6>/2) + Aa{k' 2 z + 2 ik'y/3 — z(3 2 ) sin 2 (0/2) 


( 26 ) 
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where the notations v = 2 sinh(ah); w = 2 cosh (ab);y = 2 cosh(/36);z = 2 skill (/%) have 
been used. Thus we find all the left and right handed scattering amplitudes as, 


M 2 i 

M 22 


M 


r r = 


12 


M 


U 


22 


det[M] 

M 22 




1 


with det[M] 


1 ; 


(27) 


and the transmission and reflection coefficients are obtained from these amplitudes as, 


Ti =| ti |~= T r ; Ri =| n | 2 ; Ri —\ ri 


(28) 


The phase difference between incident and transmitted waves is calculated from the trans¬ 
mitted amplitude in Eq. [271 With this phase difference the tunneling time for a single 
complex barrier is calculated by using Eq. U As a check we chose V c — 0 for which 
we should recover the tunneling time for single real barrier. At this limiting case M 22 is 
reduced (as 6 = 0 in Eq. (TT 2 D ) to, 


M 22 — 


giby/E 


2 VE^V r 


E cosh (by/V r - E ) - i(2E - V r ) sinli (by/V r - E) 

2 Ve^v~^e 


(29) 


due to which the transmission amplitudes in Eq. (127)1 is reduced to Eq. (E]) and so the 
tunneling time for a single real barrier can be re-obtained and the two channel study for 
the complex barrier is being verified. The tunneling time for a single complex barrier 
shows HF effect in the case of weak coupling (i.e. small V c ) which is already discussed in 

m ■ 

We are now at the position to calculate the tunneling time for an array [TO) SI) 02] 
of complex barriers. For that we consider an array consisting n complex barriers, each of 
width b and consecutively gaped by a length L. Each barrier potential in the array has 
been expressed in the same way as written in Eq. [7] but now the span of the (n + l) th 
potential will be decided by, 


a 1 = n(b + L) ; a 2 = b + aq ; n — 0,1, 2, 3.... (30) 

We derive the elements of individual transfer matrices for the nth barrier potential using 
Eqs. (123112(11) and by replacing the corresponding cq, a 2 as per in Eq. (l30l) . We denote the 
M-matrix for the first barrier (i.e. n = 0, span is between 0 to b) as Mq for the second 
barrier (i.e. n — 1, span is between b+L to 2b+L) as Mi and so on. In this way the total 
transfer matrix for an array of n barriers is written as, 

M tot = M n _i M 3 .M 2 .Mi.M 0 = n £=°iMj (31) 

Therefore the transmission amplitude for the array of n barriers is t tot = jjtzr- Now this 

transmission amplitude can be written also as t tot —\ t tot \ e 16 ^ where 5(k) is the phase 
difference between incident and transmitted waves for the array of barriers. Now the total 
span of the interacting potential is the sum of all the individual barriers of width b and 
their consecutive separations of length L. We can think the array of barriers equivalently 
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as a single potential distribution [HJ of width nb + (n — 1 )L. Then using Eq. fl3]) we 
obtain the tunneling time for the array of barriers as, 


dS(E ) nb + (n — 1)L 
dE hk/m 


(32) 


As the matrix multiplication in Eq. (l3Tj) is too lengthy to handle analytically we numer¬ 
ically extract the behavior of tunneling time to see the presence of HF effect in various 
situations. In the following sections we have demonstrated the results related to HF effect 
in tunneling time for array of complex barriers. 


4 Results and Discussions 

In this section we consider the expression of M tot in Eq. (T3Th to find the effect of various 
parameters on the tunneling time. We consider both absorptive as well as emissive (A is 
negative) inelastic channels in our discussion to observe HF effect with respect to number 
of barriers by varying the coupling, width of the individual barrier and separation between 
the adjacent barriers. We further discuss the tunneling time resonances in the array of 
complex barriers. For the sake of more realistic system we consider the coupling V c in a 
random manner for different barriers in one array and in another one we chose random 
excitation for different barriers. We are also able to capture the HF effect even in such 
realistic systems. 

4.1 Effect of coupling 

In evanescent mode of tunneling of waves trough an array of real barrier potentials the 
tunneling time always saturates with the number of barriers. However this is not the case 
for the complex barriers due to the presence of absorption or emission in the potential. In 
the two channel study the strength of the absorption or emission depends on the potential 
V c which couples elastic and inelastic channels. We numerically calculate the tunneling 
time by considering an array of complex barriers with same height, width and with equal 
spacing. The effect of V c on the tunneling time and on absorptivity (1 — R — T) is plotted 
in Fig. [2] We observe the saturation of tunneling time with respect to the number of 
barriers when V c is small and HF effect disappears for higher value of V c . No saturation 
of tunneling time with respect to number of barriers occurs when absorption is more in 
the array of complex barriers. Fig. [2] further shows that absorptivity also saturates with 
respect to the number of barriers. 
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tunneling time 
- absorption 
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no. of barriers 



no. of barriers 


(c) 


V =1.0 

c 



V =1.9 

c 



(b) 


(d) 


Figure 2: Tunneling times and absorptivity are shown with the number of barriers for four 
different arrays with different values of V c with E = 0.5, V r — 1, V* — 1, A = 0.02, b = 
0.2, L = 0.01. In Fig.l(a)-(c) absorption saturates almost at low value and we observe HF 
effect with respect to number of barriers. Fig. 1(d) show no HF effect for high V c (= 1.9) 
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V =0.3, b=1.0 

c 



Figure 3: Tunneling times are plotted with the number of barriers for six different arrays 
where the width of the potentials in each array is increasing with a fixed complex coupling 
V c = 0.3 and with the same values of the other parameters as in Fig. [H 

4.2 Dependency of tunneling time on width 

For a single complex barrier with weak coupling tunneling time saturates as the width of 
the barrier increases. In case of the array of the complex barrier with weak coupling we 
observe HF effect with respect to number of barrier only when width of the individual 
barrier is kept between some specific values. This has been demonstrated in Fig. [3j For 
a fixed weak value of coupling V c = 0.3 we plot the tunneling time with respect to the 
number of the barriers for different values of width b. We observe HF effect occurs for this 
particular value of V c only when width b is kept between approximately 0.1 to 0.6. This 
behavior of HF effect clearly distinct from that of array of real barrier and is demonstrated 
in Fig. |3j 

In the case of complex barriers saturation of tunneling time occurs for certain ranges 
of the width of the individual barriers. We observe HF effect for 0.2 < b < 0.6 for V c = 0.3 
keeping other parameters fixed and tunneling time starts saturating again when b > 8 
[Fig. SJa)] 0. This result is clearly distinct from the array of real barriers where tunneling 
time is independent of b except at certain small values [Fig. |4fb)] . 

4 However for exceptionally thick barriers the transmitting wave packet may be distorted. 
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V =0.0 
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Figure 4: In (a) and (b) the tunneling times with the width of barriers are plotted for an 
array of complex barriers ( V c = 0.3) and an array of real barriers ( V c = 0) respectively by 
keeping all other parameters same as in Fig. 0 The insets of\^fa) zoom out the behavior 
of tunneling time for smaller values of width (red line) and lager values of width (green 
line). The same range of width has been enlarged in the inset of\J^b) to capture distinct 
behavior of tunneling time for real array. 
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4.3 Tunneling time resonances 


The resonances in tunneling time have been noticed for the case of real barriers [37] . In 
case of one or two real barrier HP it has been reported that for certain values of the 
separation between adjacent barriers the wave packet never emerges indicating a very 
high value in tunneling time. In this subsection we numerically study such resonances 
in the case of array of complex barriers. In particular we show how resonances in tun¬ 
neling time changes with coupling between elastic and inelastic channels and separation 
between adjacent barriers. We observe that resonances in tunneling time for evanescent 
wave gradually disappear if we increase V c coupling between the elastic and inelastic 
channels. This is explained clearly in Fig. [5] by changing V c in the array of total 20 real 
barriers (Fig. [5](a)) and 20 complex barriers (Fig. [5](b, c, d)). To see the effect of complex 
coupling on the resonances we chose a much opaque barrier by increasing both real and 
imaginary parts of the barrier height (i.e V r , Vj) up to a value so that we can increase 
the value of V c without disturbing the evanescent mode of tunneling. The first resonance 
disappears when complex coupling V c raise from 1 to 14. Unlike the case of real barriers 
the resonances in the case of array of complex barrier is regulated by increasing inelasticity. 
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Figure 5: Tunneling times are shown with consecutive separation of barriers in the array. 
Figure (a) is for an array of real barriers and (b,c,d) are for arrays of complex barriers 
with different complex coupling V c = 0, V c = 5, V c = 8 and V c = 14 respectively with the 
other parameters same as in Fig. [H Resonances is regulated as we increase the inelasticity 
in the system. 
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4.4 HF effect in randomly arranged arrays 


In this section we would like to point out HF effect is not an artifact of regularity in 
terms of coupling V c and/ or A. In a realistic system the inelasticity may be different in 
different barriers. Keeping this in view we chose the coupling V c randomly for different 
complex barriers in the array and study tunneling time by considering arbitrary excited 
state energies (A). We show HF effect exists even for such systems and are demonstrated 
in Fig. [6] (a, b). Similarly we consider another array of barriers in which individual 
barriers correspond to different inelastic channel with random values of A. We show the 
saturation in tunneling time in such a array of barrier for the propagation of wave packet 
in the evanescent mode. This is graphically represented in Fig. O (c, d) with keeping the 
other parameters same as Fig. [2l 


14 


0.35 

( 

0.3 

0.25 

>° 0.2 
0.15 9 


0.1 

0.05 


°°o o cP °o? oo °0 0 (9° 

° o ® 0 °o° ° O c 

°%°o °o° ° 


o o 

9? O o °o ° ° ° 

cr _ o c 


o co 


9c# ° 


° 9) 


60 80 100 120 140 160 

no. of barriers 


~o~ 


o 9> 


cb 


o 


o 


O O o ° ° 

( 

o c> 


cP 


o 


o 


O o 

I 

o O 


(9 o O 
° O 00 O O 

° 


cb - 


oo 


o W O ° ° ° 

° <§>§ n ° O ° ° 

8 ° O O oo o o 
O o ^ ° Q o 
o ° ° O o o 


Cb 


—O— 


o 


o 


o 

o o 


o o 


60 80 100 120 140 160 

no. of barriers 


0.01 > V c <0.4 ;A=0.2 




(b) 


(d) 


Figure 6: The random distributions (between a range of 0.01 to O.f) of coupling V c and 
excitation energy A of each barriers in two different arrays have been shown in (a) and 
(c). The corresponding tunneling times for these two arrays are shown in the figures (b) 
(for fixed A = 0.02 ) and (d) (for fixed V c = 0.3 ) respectively. 


4.5 Emissive inelastic channel 

In the previous subsections we consider absorptive (A is positive) inelastic channel to 
demonstrate HF effect. Now we concentrate emissive inelastic channel by considering 
target in an excited state A. The Schroedinger equation in this inelastic channel is then 
written by replacing A by —A in Eq. (JH]). Therefore an emission is being associated 
with the interaction of waves and the potential. For E « A both elastic and inelastic 
processes occur in presence of such a target. We observe HF effect with respect to number 
of barriers at particular very low energy of the incident wave packet (see Fig. [TJ) . In the 
inset of each graphs of Fig. [7] we have plotted the relative tunneling time ( T10 °~ r ) with 
respect to an array of 100 barriers. Fig. [7](b) shows the HF effect at around E = .00035 
for A = 3.6. Occurrence of HF effect at particular low incident energies is an interesting 
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Figure 7: Tunneling times with the number of barriers in an array for fixed A = —3.6 
(with V r = 10, V c = 14 and other values are same as Fig. d|) are shown when the incident 
particle energy is considered very small, from a value of E = 0.00001 to E = .0008. 

result in the case of array of complex barriers and need further investigation. 


5 Conclusions 

We have studied the HF effect for an array of complex barrier potentials to unfold various 
new features associated with this interesting effect by using stationary phase method. We 
have constructed the total scattering transfer matrix by multiplying the transfer matrices 
for the individual barriers to calculate the tunneling time for a wave packet through such 
an array of barriers. The tunneling time saturates with respect to the number of barriers 
depending on the different parametric values in the system. Saturation crucially depends 
on the coupling potential V c which couples elastic and inelastic channels of propagation 
when other parameters are held fixed. We have observed HF effect only for low absorption 
(i.e. small V c ) in the system. Saturation in tunneling time with respect to number of 
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barriers are observed for certain ranges of the width of barrier unlike the situation for the 
real array where saturation is always achieved beyond a certain value of the width. In 
case of real array for certain values of the width of the barrier and separation of adjacent 
barrier the tunneling time is extremely high and the wave packet never emerges from such 
barriers. We have shown such a resonance behavior in tunneling time is regulated in case 
of array of complex barriers. For the sake of realistic systems we further have studied 
the HF effect for the array of barriers with random value of coupling with arbitrary value 
of excited energy (A) and with fixed inelasticity with random (A). In both cases we 
have shown the saturation of tunneling time with respect to number of barriers. Finally 
we have observed HF effect even for the case of emissive inelastic channel. Surprisingly 
saturation of tunneling time occurs there at some particular very low incident energy. 
This needs further investigations. 
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